Interface solitons in quadratically nonlinear photonic lattices 
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We study the properties of two-color nonlinear localized modes which may exist at the interfaces 
separating two different periodic photonic lattices in quadratic media, focussing on the impact of 
phase mismatch of the photonic lattices on the properties, stability, and threshold power require- 
ments for the generation of interface localized modes. We employ both an effective discrete model 
and continuum model with periodic potential and find good qualitative agreement between both 
models. Dynamics excitation of interface modes shows that, a two-color interface twisted mode 
splits into two beams with different escaping angles and carrying different energies when entering a 
uniform medium from the quadratic photonic lattice. The output position and energy contents of 
each two-color interface solitons can be controlled by judicious tuning of the lattice parameters. 

PACS numbers: 42.65.-k, 42.65.Tg, 42.65.Wi 



I. INTRODUCTION 

Electromagnetic surface waves are waves localized at 
an interface separating either two homogeneous (one of 
them has to be surface-active, i.e., exhibiting a negative 
permittivity [1]) or homogeneous and periodic dielectric 
media [2]. In addition, nonlinear dielectric media can 
support nonlinear guided waves localized at or near the 
surfaces, and different types of nonlinear guided waves in 
planar waveguides have been studied extensively about 
20 years ago [21 H] • Recently, the interest in the study of 
electromagnetic surface waves has been renewed after the 
theoretical prediction [5] and subsequent experimental 
demonstration [6] of nonlinearity-induced self-trapping 
of light near the edge of a one-dimensional waveguide 
array with self-focusing nonlinearity, that can lead to 
the formation of a discrete surface soliton. A related 
effect of light localization and the formation of surface 
gap solitons have been predicted theoretically and ob- 
served experimentally for defocusing periodic nonlinear 
media 0IE]. In addition, the concept of nonlinear sur- 
face and gap solitons has been extended to the case of 
an interface separating two different nonlinear periodic 
media [SB ED E2] • 

Surface solitons are usually considered for one- 
frequency modes propagating in cubic or saturable non- 
linear media. However, multicolor discrete solitons in 
quadratically nonlinear lattices have been studied theo- 
retically in both one- and two-dimensional lattices |13l 
ITU [T5J [TB] irrespective to the surface localization effects. 
Only Siviloglou et al. [T7j studied discrete quadratic sur- 
face solitons experimentally in periodically poled lithium 
niobate waveguide arrays, and they employed a discrete 
model with decoupled waveguides at the second harmon- 
ics to model some of the effects observed experimentally. 

More elaborated theory of one-dimensional surface 
solitons in truncated quadratically nonlinear photonic 
lattices, the so-called two-color surface lattice solitons, 
has been developed recently by Xu and Kivshar [18] 
who analyzed the impact of the phase mismatch on the 



existence and stability of nonlinear parametrically cou- 
pled surface modes, and also found novel classes of one- 
dimensional two-color twisted surface solitons which are 
stable in a large domain of their existence. 

The purpose of this paper is twofold. First, we extend 
the analysis of two-color surface solitons to the case of two 
semi-infinite photonic lattices with quadratic nonlineari- 
ties. We study, for the first time to our knowledge, two- 
color interface solitons in photonic lattices with quadratic 
nonlinear response. We analyze the effect of mismatch 
on the existence, stability, and generation of such novel 
surface states. Second, for the analysis outlined above 
we employ two different approaches widely used in the 
literature: The coupled-mode theory described by the 
discrete parametrically coupled equations for the funda- 
mental and second-harmonic fields, and also the contin- 
uous model with a periodic potential, and demonstrate 
that both models give the same qualitative results. 

The paper is organized as follows. In Sec. II we discuss 
the two-color interface localized modes in the framework 
of a discrete model where the interface is modeled by 
a jump of the propagation constant, assuming that the 
matching conditions are nearly satisfied for both semi- 
infinite lattices. Section III is devoted to the analysis 
of a more general case described by a continuous model 
with a periodically varying potential. Thus we consider 
the light beam propagating along the interface between 
two dissimilar optical lattices imprinted in quadratic non- 
linear media. Besides the study of the properties of inter- 
face solitons, we demonstrate the manipulation of surface 
solitons by tuning the lattice and waveguide parameters. 
Finally, Sec. IV concludes the paper. 



II. DISCRETE MODEL 

We consider the propagation of light in a one- 
dimensional photonic lattice of a finite extent imprinted 
in quadratic nonlinear media, which involves the inter- 
action between fundamental frequency (FF) and second- 
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harmonic (SH) waves. Light propagation is described by 
the following coupled nonlinear discrete equations [131118] 

du 

i—^ + C u (u n+1 + u n -i) + 2ju* n v n exp(+if3z) = 0, 

.dvn + c v (v n+ i + u n _i) + ju^ exp(-i/3z) = 0, (1) 
az 

where u n and v n are the normalized amplitudes of the 
FF and SH waves, respectively, C u and C v are the cou- 
pling coefficients, 7 characterizes the second-order non- 
linearity, and (3 is the effective mismatch between the two 
harmonics. 

We look for stationary two-mode solutions of Eq. in 
the form, u n (z) = U n exp(ibz) and v n {z) = V n exp(2ibz — 
i/3z), where b is the propagation constant, and obtain the 
nonlinear algebraic equations for the (real) mode ampli- 
tudes, U n and V„, 

-bU n + C u (U n+1 + U n - 1 ) + 2 1 U n V n = 0, 
-2bV n + C v {V n+ x + V n -x) + p n V n + 1 Ul = 0, (2) 

where j3 n — /3i on the left side of the interface, and 
ftn = 02 ° n the right side. We solve Eqs. Q nu- 
merically, via the well-known multi-dimensional Newton- 
Raphson method, using the results of the anti-continuum 
limit as initial conditions for the algorithm: When b 3> 
Max[/3i;/3a], it is easy to obtain, V n « b/{2j) and 
U n 6/7. Figures l(a,b) show an example of the sim- 
plest interface mode, where the fields propagating along 
the waveguides are presented as a superposition of the 
waveguide modes, U(x) = ^2 n U n <j)^^(x — n), where 
<f>(ff}(x) is the single FF mode centered on guide n; sim- 
ilarly for V{x). Using the same numerical approach and 
starting from the anti-continuum limit, we find several 
other families of two-color localized modes located at the 
interface. These modes provide a generalization of the 
surface modes known for truncated one-dimensional lat- 
tices located at different distances from the edge, and cor- 
responding to a crossover between the interface and bulk 
discrete solitons as discussed earlier [T^]. In addition, we 
find novel classes of the so-called interface twisted modes, 
and an example of such mode is shown in Figs. [TJc,d). In 
general the spatial profile of these modes are asymmetri- 
cal for /3i 7^ /?a. 

For given values of C u , C v , (3\, and /3a, we compute 
the minimum value of the propagation constant b and 
the power 



P = P u + P v = Y,(K\ 2 + 2\v n \ 2 ), 



(3) 



needed to create a localized mode at one of the inter- 
face sites. Figures 2(a,b) show examples of P m ; n curves 
for two cases. In the first case, when no interface is 
present, there is no minimum power to generate a lo- 
calized mode for small values of (3. Outside this interval, 
power increases more or less symmetrically with increase 
in \ f3\. On the other hand, in the presence of an interface 



we observe that the generation of an interface localized 
mode require a finite power, as shown in Fig. 2(b). The 
numerically-observed stability of these modes seems to be 
in agreement with the stability results of the continuous 
model (see section III). 
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FIG. 1: Examples of two-color interface solitons for C« = 1, 
C v = 0.5 and (a,b) 6 = 5, Pi = 3, = -3, and (c,d) 6 = 4, 
Pi = 3, @2 = —3. Shown are the amplitudes of the fundamen- 
tal (a,c) and second-harmonic (b,d) components. Case (a,b) 
corresponds to a fundamental interface mode, while (c,d) cor- 
responds to an interface "twisted" mode. 
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FIG. 2: Minimum power P m i n of the simplest two-color inter- 
face soliton for C u = 1, C v = 0.5, and (a) Pi = Pi = P, (b) 
Pi = and p 2 = p. 

The dynamical excitation of the interface mode can be 
achieved by launching all of the initial power in one of the 
guides adjacent to the interface. When both, FF and SH 
fields are initially present, a localized mode is formed if 
the power is strong enough. In the more experimentally 
relevant case when only the FF field is initially present, 
the threshold power increases significantly. After the in- 
terface mode has been established, the power exchange 
between the FF and SH fields proceeds along the longitu- 
dinal direction with a spatial period that decreases with 
increasing initial power, for Pi, 02 fixed. On the other 
hand, for a fixed value of initial power in the FF field, 
enough to excite the interface mode, the spatial period 
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for power exchange between the fields decreases with an 
increase in the value of at the initial launching guide. 
Figure [3] shows an example of the dynamical excitation 
of a simple odd mode and a twisted mode, where only 
the FF field is excited initially. Internal oscillations of 
the field as well as asymmetrical power sharing for the 
twisted mode are apparent. 
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FIG. 3: Dynamical excitation of interface modes. Power 
content in the FF field for (a) Simple odd mode and (b) 
Twisted mode. Only FF field is excited initially. Here 
C u = 1, C v = 0.5, pi = 3, /3 2 = -3. 



III. CONTINUOUS MODEL WITH A 
PERIODIC POTENTIAL 

In this section, we extend the analysis to a more 
general model (a continuous model with a periodically 
varying potential) to describe propagation of light at 
the interface between two dissimilar optical lattices im- 
printed in quadratic nonlinear media, which involves the 
interaction between fundamental frequency and second- 
harmonic waves. Light propagation is described by the 
following coupled nonlinear equations [18 



.du d\ d 2 u 



' dz 
,dv 
dz 



2 dx 2 
Ydx 2 



— u*vexp(—if3z) — pR(x)u, 

— u 2 exp(i(3z) — 2pR(x)v, 



(4) 



where u and v represent the normalized complex am- 
plitudes of the FF and SH fields, x and z stand for 
the normalized transverse and longitudinal coordinates, 
respectively, (3 is the phase mismatch, and d\ = —1, 
d 2 = —0.5; p is the lattice depth; the function R(x) — 
s[l — cos(K 2 x)] at x < and R(x) = 1 — cos(A"!2;) 
at x > describes the profile of a waveguide array 
formed with two dissimilar lattices with modulation fre- 
quencies K\ and K 2 , respectively, where s is the relative 
lattice depth for the left-side one [a typical profile for 
such tunable waveguide array is shown in Fig. [2]. In 
typical quadratic nonlinear crystals, for a beam width 
of ~ 15/im, the distance z in the range — 30 corre- 
sponds to a few centimeters, and the peak intensities will 
be in the range of 0.1 — lOGW/cm 2 for the formation of 



lattice solitons at wavelengths A = Ifm; a refractive in- 
dex modulation depth of the order of 10 -4 corresponds 
to the lattice depth p ~ 1 [18 . The system of Eq. Q 
admits several conserved quantities including the power 
P = C OQ {\u\ 2 + \v\ 2 )dx. 
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FIG. 4: Interface structure created by two different quadratic 
lattices with tunable guiding parameters. 

The stationary solutions for the lattice-supported in- 
terface solitons can be found in the form u(x, z) — 
Wi(x)exp(ib\z) and v(x, z) — w 2 (x)exp(ib 2 z), where 
^1,2(2^) are real functions, and b\^ 2 are real propagation 
constants satisfying b 2 — (3 + 2bi . Families of interface 
solitons are determined by the propagation constant b±, 
the lattice depth p and s, and the phase mismatch (3. For 
simplicity, we set the modulation parameter K\ = 4 and 
vary K 2 . To analyze stability we examine perturbed solu- 
tions u{x,z) — [wi(x) + Ui(x, z) + iVi(x 1 z)]exp(ibiz) and 
v(x, z) — [w 2 {x) + U 2 (x,z) + iV 2 {x, z)]exTp(ib 2 z), where 
real parts Ui, 2 and imaginary parts Vi. 2 of perturba- 
tion can grow with complex rate S. The linearization 
of Eq. Q around w±. 2 yields the following eigenvalue 
problem: 



8U1 = 


di d 2 V x 
2 dx 2 


{wiV 2 - 


w 2 V 1 )-pRV 1 + b 1 V 1 , 


5V 1 = 


d 1 d 2 U x 
2 dx 2 


+ {w x U 2 


+ w 2 U 1 )+pRU 1 -biUx, 


6U 2 = 


d 2 d 2 v 2 

2 dx 2 


2w\Vi - 


2pRV 2 + b 2 V 2 , 


5V 2 = 


d 2 d 2 u 2 

2 dx 2 


+ 2wiUi 


+ 2 P RU 2 -b 2 U 2 , (5) 



which we solve numerically to find the growth rate 6. 

The tunable waveguide array shown in Fig. [4] supports 
rich families of two-color interface lattice solitons. In 
semi-infinite waveguide arrays (s = 0) the properties 
of two-color surface solitons have been investigated [TB] • 
Here we are interested in the case s ^ thus we can show 
how to manipulate two-color lattice interface solitons by 
tuning the lattice and waveguide parameters s and modu- 
lation frequency K 2 . The simplest example of a two-color 
interface mode (odd soliton) is shown in Fig. [5]ja), from 
which one notices that the peak of this mode coincides 
with one of the maxima of the lattices. As shown in 
Fig. [5^b) its total power is almost a monotonic function 
of the propagation constant. It it important to note that 
the existence domain of two-color interface solitons is de- 
termined by the waveguide parameters, thus the cutoff 
value of propagation constant varies depending on the 
lattice period [Fig. [5][c)] and depth [Fig. f5jd)]. A lin- 
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FIG. 5: (Color online) Two-color odd interface solitons (a) 
FF (black solid curve) and SH (red dashed curve) field dis- 
tributions of two-color odd interface solitons with bi = 4 at 
s = 0.5 and K2 = 2.6. (b) Power versus propagation constant 
for different parameters as marked in the plot, In (a,b) phase 
matching /3 = 0. (c, d) Cutoffs of the propagation constant 61 
versus (c) the lattice period K2 at s — 0.5 and (d) the relative 
lattice depth s at K% = 3 for different /3. 
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FIG. 6: (Color online) Two-color twisted interface solitons 
(a) FF (black solid curve) and SH (red dashed curve) field 
distributions with 61 = 4 at s — 0.5 and K2 = 2.6. (b) 
Power versus propagation constant for different parameters 
as marked in the plot [Solid and dotted curves correspond to 
stable and unstable branches, respectively], (c, d) Cutoffs of 
the propagation constant b\ versus (c) the lattice period K2 
at s — 0.5 and (d) the relative lattice depth s at K2 = 3. Here 
phase matching /3 = 0. 



ear stability analysis reveals that these odd solitons are 
stable almost in the entire domain of their existence. 

Besides the odd solitons described above, we find that 
the tunable waveguide with two dissimilar optical lat- 
tices also supports higher-order interface lattice solitons, 
which can be viewed as combinations of several odd soli- 
tons in engineered phases. Fig.[6ja) shows the case where 
the FF field features the out-of-phase combination be- 
cause the in-phase combinations are unstable; we term 
the mode with an out-of-phase combination as "twisted 
interface lattice solitons" . There exist cutoff values of 
propagation constant for twisted interface solitons, which 
depend on the lattice period [Fig. [6|c)] and lattice depth 
[Fig. (6^d)]. Note that the profiles of twisted modes are 
asymmetric due to the dissimilarity of the lattices on 
two sides [Fig. [6^ a)]. An important result is that such 
asymmetric twisted modes are stable when the power ex- 
ceeds a critical value [unstable region is shown in dotted 
curves in Fig. |6jb)], which is confirmed by both, our lin- 
ear stability analysis and the direct numerical simulation 
of the soliton generation from the fundamental beam in 
the framework of the model Q . Generation and stable 
propagation of asymmetric twisted surface solitons is also 
shown in Fig. (7^b). 

One of the key issues in this paper is to demonstrate 
the control and manipulation of two-color interface soli- 
tons by tuning the guiding parameters of the waveguide 
array. Removal of the lattices causes each component 
of the twisted soliton to fly apart because the inter- 



action force between neighboring constituents is repul- 
sive, due to the out-of-phase combination that forms a 
twisted mode. Figs. |7jc,d) show a scenario in which a 
two-color interface twisted soliton splits into two beams 
carrying different energies when entering into a uniform 
medium from the quadratic photonic lattices. Note that 
the power carried by each beam is different due to the 
lattice-induced asymmetric profiles of two constituents of 
the twisted solitons and as a result, the escaping angles 
are different for the two beams. Thus, we can control the 
output position of the two beams and the power sharing 
by tuning the guiding parameters of the quadratic pho- 
tonic lattices. The whole manipulation and switching 
scenarios of twisted two-color interface solitons in tun- 
able lattices are summarized in Fig. [8] from which one 
can see clearly that the escape angle and energies of two 
constituents of twisted solitons can be controlled by tun- 
ing the relative lattice depth s [Figs.[8la,b)] and the mod- 
ulation frequency K2 of the lattices [Figs. [8^c,d)]. These 
results suggest novel opportunities to control the flow of 
light in tunable photonic systems. 



IV. CONCLUSIONS 

We have studied parametric localization of light at an 
interface separating two different quadratically nonlin- 
ear photonic lattices and determined the conditions for 



5 




64 



CD 
O 

o3 32 

CO 

s 



V 


/ 


] 





64 



-10 -8 8 
X 




the existence of interface localized modes. We have ana- 
lyzed the impact of the phase mismatch on the properties 
and stability of interface localized modes, as well as the 
threshold power for their generation. We have employed 
two different approaches: the coupled-mode theory and a 
continuous model with a periodic potential. Both models 
give basically the same qualitative results. Our results 
reveal that asymmetric phase-twisted two-color surface 
solitons are stable in a broad region of system parame- 
ters, but they can split into two solitons when entering a 
homogeneous medium. We have demonstrated that the 
output position and energies of constituents of the two- 
color interface solitons can be controlled by tuning the 
lattice and waveguide parameters. 
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FIG. 7: Excitation of interface solitons from only FF field 
as input beam: FF field distributions of generated odd (a) 
and twisted (b) interface solitons. (c,d) Splitting of generated 
twisted interface solitons at the boundary between a uniform 
and periodic media for (c) FF and (d) SH fields. Here = 0, 
s = 0.5, and Ki=K 2 = 4. 
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FIG. 8: (a) Escape angles and (b) energy sharing of interface 
twisted solitons as a function of parameter s for K2 = 4. 
(c) Escape angles and (d) energy sharing of interface twisted 
solitons as a function of parameter K2 at s = 0.5. Here 
Phase matching f3 = 0. In all cases solid and dashed curves 
correspond to the left and right beams in the twisted soliton, 
respectively. 
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